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neutrinos produced by photo-meson (or p-p) interactions in sources of 
size not much larger than the proton photo-meson (or p-p) mean-free- 
path. This bound applies, in particular, to neutrino production by either 
AGN jets or GRBs. The upper limit is two orders of magnitude below 
the intensity predicted in some popular AGN jet models and therefore 
CN| ■ contradicts the theory that the cosmic gamma-ray background is due 

to photo-pion interactions in AGN jets. The upper bound is consistent 
qq ■ with our predictions from GRB models. The predicted intensity from 

GRBs is E 2 dN/dE ~ 0.3 x 10^ 8 GeV/cm 2 ssr for 10 14 eV < E < 10 16 eV; 
^ ; we also derive the expected intensity at higher energy. 

. PACS numbers: 95.85.Ry, 98.70.Rz, 98.70.Sa, 14.60.Pq 



<P . I. INTRODUCTION 



Large area, high-energy neutrino telescopes are being constructed to detect cosmologi- 
cally distant neutrino sources (see [fJJ for a review). The main motivation of a search for 
cosmological high-energy neutrino sources derives from the fact that the cosmic-ray energy 
spectrum extends to > 10 20 eV and is most likely dominated above ~ 3 x 10 18 eV by an extra- 
Galactic source of protons @. High-energy neutrino production is likely to be associated 
with the production of high-energy protons, through the decay of charged pions produced 
by photo-meson interaction of the high-energy protons with the radiation field of the source. 
Gamma-ray bursts (GRB) || and active galactic nuclei (AGN) jets ||] have been suggested 
as possible sources of high-energy neutrinos that are associated with high-energy cosmic- 
rays. The predicted neutrino fluxes may be detectable with high-energy neutrino telescopes 
of effective area ~ 1 km 2 . 

We show here that high-energy cosmic-ray observations set a model-independent upper 
bound to the expected high-energy neutrino fluxes and are in conflict with the theory that 
the cosmic gamma-ray background is due to photo-pion interactions in AGN jets. The 
upper bound is stated in Eq. (Rf) and is illustrated in Fig. 1. The demonstration that the 
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AGN jet models for the gamma-ray background are in conflict with the cosmic ray data is 
given in Section III. We also give a more detailed prediction than we gave in ref. || for the 
expected GRB neutrino spectrum and discuss the compatibility of the detailed results with 
the general secondary-particle cooling constraints derived by Rachen & Meszaros Q. 

It has been suggested that neutrinos may be produced in the cores of AGNs (rather than 
in the jets) by photo-meson interaction of protons accelerated to high energy in the AGN 
core; in this scenario neutrinos are produced very close to the central black-hole ||. In 
this model, the proton photo-meson optical depth is very high, r P7 ~ 100, and high-energy 
nucleons do not escape the source. By construction, there can be no observational evidence 
except neutrinos for, or against, the hypothesized luminous AGN accelerator of high-energy 
protons. The hypothesized black-hole accelerators are "neutrino-only" factories. Therefore, 
cosmic-ray observations can not set a limit to neutrino emission in this model. On the 
other-hand, this model can not explain, and is therefore not supported by, the existence of 
the extra-Galactic high-energy cosmic-ray flux. 

This paper is organized as follows. In section II we derive the general upper bound to 
neutrino fluxes from p — 7 interactions for sources optically thin to p — 7 reactions. We 
compare in this section the upper limit to the predictions from different models for neutrino 
sources. We also show that the upper bound cannot be avoided by cosmological evolutionary 
effects (§IIc) or by invoking magnetic fields (§111). In section IV we discuss the implications 
of the upper bound for AGN jet models of neutrino production. In section V we discuss the 
implications for the GRB model of high energy neutrino production, derive a more detailed 
prediction (compared to our prediction in ref. ||) for the expected GRB neutrino spectrum, 
and compare our results with those of other authors. We discuss in Section VI our main 
conclusions. 

II. UPPER BOUND TO THE NEUTRINO FLUX 

We first derive in section IIA the upper bound to the high-energy neutrino flux from the 
sources at redshift z < 1 that produce the observed cosmic-rays at energies greater than 
10 18 eV. We compare in section IIB the upper bound with the predictions of current models. 
In section IIC we discuss the modification of the upper bound by unobserved sources of 
cosmic rays at larger redshift. 

A. Derivation of the upper bound 

Cosmic-ray observations above 10 17 eV indicate that an extra-Galactic source of protons 
dominates the cosmic-ray flux above ~ 3 x 10 18 eV ||, while the flux at lower energies is 
dominated by heavy ions, most likely of Galactic origin. The observed energy spectrum 
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of the extra-Galactic component is consistent with that expected for a cosmological dis- 
tribution of sources of protons, with injection spectrum dN C R/dE C R oc Eq r , as typically 
expected for Fermi acceleration M. The energy production rate of protons in the energy 
range 10 19 eV to 10 21 eV is Eq R ' 10 ' ~ 5 x 10 44 erg Mpc _3 yr _1 M, if the observed flux of 
ultra-high-energy cosmic-rays is produced by sources that are cosmologically distributed. 
The energy-dependent generation rate of cosmic-rays is therefore given by 

J AT .[10 19 ,10 21 ] 

= Rfe - ^V g Mpc-V'. (!) 

If the high-energy protons produced by the extra-galactic sources lose a fraction e < 1 of their 
energy through photo-meson production of pions before escaping the source, the resulting 
present-day energy density of muon neutrinos is E^dN^/ ' dE v « 0.25etHEc R dNcR/dEcR, 
where tjj ~ 10 10 yr is the Hubble time. For energy independent e the neutrino spectrum 
follows the proton generation spectrum, since the fraction of the proton energy carried by a 
neutrino produced through a photo-meson interaction, E v ~ 0.05E P , is independent of the 
proton energy. The 0.25 factor arises because neutral pions, which do not produce neutrinos, 
are produced with roughly equal probability with charged pions, and because in the decay 
7r + — ► fi + + Ufj, —>■ e + + v e + Up + Up muon neutrinos carry approximately half the charged 
pion energy. Defining 7 max as the muon neutrino intensity [y^ and combined) obtained 
for e = 1, 

J max w 0.25Z z t H ^-E 2 CR ^^ w 1.5 x KT^zGeVcm-V^r- 1 , (2) 
47r dE C R 

the expected neutrino intensities are 

o c ndN v 1 . , 

El^ = -^K^ = ^el max , (3) 

The quantity (,z in Eq. (|) is of order unity and has been introduce here to describe the 
possible contribution of so far unobserved high redshift sources of high-energy cosmic rays 
and to include the effect of the redshift in neutrino energy. We estimate in section IIC 



B. Upper bound versus current models 

Figure 1 compares the neutrino intensity predictions of GRB and AGN jet models with 
the intensity given by Eq. (Q). The AGN core model predictions are shown for completeness. 
The intensities predicted by both AGN jet and core models exceed J max by typically two 
orders of magnitude. 
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FIG. 1. Comparison of muon neutrino intensities (i>u and P„ combined) predicted by different 
models with the upper bound implied by cosmic ray observations. The dash-dot lines give the 
upper bound Eq. (|2[) corrected for neutrino energy loss due to redshift and for possible redshift 
evolution of the cosmic-ray generation rate. The lower line is obtained assuming no evolution, and 
the upper line assuming rapid evolution similar to the evolution of the QSO luminosity density. The 
AGN jet model predictions are taken from ref. |Q] (labeled 'Jetl' and 'Jet2'). The GRB intensity 
is based on the estimate presented in this paper, following ||. The AGN hidden-core conjecture, 
which produces only neutrinos and to which the upper bound does not apply, is taken from |J. 

The intensity J max is an upper bound to the intensity of high-energy neutrinos produced 
by photo-meson interaction in sources of size not much larger than the proton photo-meson 
mean-free-path. Higher neutrino intensities from such sources would imply proton fluxes 
higher than observed in cosmic-ray detectors. Clearly, higher neutrino intensities may be 
produced by sources where the proton photo-meson "optical depth" is much higher than 
unity, in which case only the neutrinos escape the source. However, the existence of such 
sources cannot be motivated by the observed high-energy cosmic-ray flux or by any observed 
electromagnetic radiation. We therefore refer in Fig. 1 to models with r 7P > 1 as "hidden 
core" models. 



C. Evolution and redshift losses 

In the derivation of Eq. (0) we have neglected the redshift energy loss of neutrinos 
produced at cosmic time t <tn, and implicitly assumed that the cosmic-ray generation rate 
per unit (comoving) volume is independent of cosmic time. The generation rate may have 
been higher at earlier times, i.e. at high redshift. Cosmic rays above 10 18 eV must originate 
from sources at z < 1. Energy loss due to redshift and pair production in interaction with 
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the microwave background implies that in order to be observed with energy E > 10 18 eV, 
a proton should have been produced at z — 1 with energy exceeding the threshold for 
photo-meson production in interaction with the microwave background at that redshift. 
Photo-meson energy loss of protons produced above the threshold would reduce the proton 
energy to the threshold value in a short time, so that their observed energy (i.e. at z = 0) 
would be ~ 10 18 eV. Thus, the cosmic ray energy generation rate given in Eq. (|1|) is the 
present (i.e low redshift, z < 1) generation rate. An increase in the cosmic-ray energy 
generation rate per unit (comoving) volume above the value of Eq. ([I]) at large redshift, 
z > 1, is consistent with observations, since it would not affect the observed flux above 
~ 10 18 eV, and the contribution from z > 1 sources to the observed flux below ~ 10 18 eV 
may be hard to detect due to the "background" of Galactic sources of heavy ion cosmic-rays 
which are most likely dominating the flux at this energy [fjf|. 

Let us consider the possible modification of Eq. (fj) due to evolution and redshift losses. 
A neutrino with observed energy E must be produced at redshift z with energy (1 + z)E. 
Thus, the present number density of neutrinos above energy E is given by 

/■Zmax dt Pmax dt 

n u {>E)= dz—h 1/ [>(l + z)E,z]=h (>E) dz—{l + z)- l f{z). (4) 
jo dz Jo dz 

Here we have used the fact that n v {> E) oc E~~ x and denoted the ratio of (comoving) 

neutrino production rate at redshift z to the present rate, no, by f{z). Comparing Eqs. 

(0) and (||), and noting that t H = J °° dz(dt/dz), we find that the intensity J max of Eq. (0) 

should be multiplied by a correction factor 

_ / 2 — dzg(z)(l + z)- 7 ^f(z) 

KZ JS°dzg(z)(l + z)-*/* ' U 

Here, g(z) = — Hq(1 + z) 5 ^ 2 (dt/dz) is a weak function of redshift and cosmology; g(z) = 1 

for a flat universe with zero cosmo logical constant. Let us assume that the neutrino energy 

generation rate evolves rapidly with redshift, following the luminosity density evolution of 

QSOs 0, which may be described as f(z) = (1 + z) a with a ~ 3 [|10j at low redshift, 



z < 1.9, f(z) = Const, for 1.9 < z < 2.7, and an exponential decay at z > 2.7 [|TTJ. Using 



this functional form of f{z), which is also similar to that describing the evolution of 
star formation rate ||12|| , we find that £z ~ 3 (with weak dependence on cosmology). For no 
evolution, f(z) = Const., we have £z ~ 0.6 < 1 (with only a weak dependence on cosmology) 
due to a redshift energy loss of neutrinos. 



III. CAN THE UPPER BOUND BE AVOIDED BY INVOKING MAGNETIC 

FIELDS? 

One may try to invent arguments in order to avoid the upper bound on the neutrino flux 
by hypothesizing strong magnetic fields, which would affect the propagation of cosmic-ray 
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protons, either in the neutrino source or in the inter-galactic medium between the source 
and the Earth. We show in this section that magnetic fields in the source can not affect 
the upper bound, and that observational constraints on inter-galactic magnetic fields imply 
that proposed field scenarios also do not affect the upper bound. 

A. Magnetic fields in the neutrino source 

One might try to argue that the upper bound J max can be avoided even for sources 
with small photo-meson optical depth, if protons are prevented from escaping the source 
by magnetic confinement. However, a photo- meson interaction producing a charged pion 
also converts the proton to a neutron, which is not magnetically confined and will escape 
a source with small photo-meson optical depth before decaying to a proton. A neutron of 
high energy E propagates a distance 100(-E/10 19 eV)kpc prior to its decay. Thus, magnetic 
fields within the neutrino source cannot be used to evade the upper bound given in Eq. (^j). 

B. Uniformly distributed inter-galactic magnetic fields 

The existence of uniformly distributed inter-galactic magnetic fields may limit the prop- 
agation distance of protons and therefore prevent their arrival at Earth from distant sources. 
However, imposing a limit on the propagation distance of cosmic-rays will not affect any of 
the arguments presented in this paper. Compare, e.g., the case where protons propagate on 
straight lines to the case where cosmic rays are confined by magnetic fields to their place of 
origin. In the former case, the present proton flux is obtained by integrating the contribu- 
tion of distant sources over redshift, while in the latter case it is obtained by integrating the 
contribution from local sources over cosmic time. For a homogeneous universe both proce- 
dures yield the same result [e.g. the integration in Eq. @ may be interpreted as integration 
over redshift or over time]. Thus, limiting the cosmic-ray propagation distance may affect 
the upper bound J max only if propagation is limited to a distance d over which the cosmic 
ray production rate, averaged over a Hubble time, is inhomogeneous. In this case, if the 
cosmic-ray production rate in our local region of size d is lower than the universe average, 
the observed cosmic-ray flux would be lower than average and the neutrino flux, which is 
homogeneous throughout the universe, may exceed a bound based on the local cosmic-ray 
flux. However, the magnetic fields required to make such a scenario viable, i.e. to confine 
cosmic rays to small enough d to significantly affect J max , are large and are inconsistent with 
observations. 

Consider a proton of energy E propagating through an intergalactic magnetic field of 
strength B and correlation length A. Propagating a distance A the proton is deflected 
by an angle ~ A/i?^, where Rl = E/eB is the Larmor radius. For the parameters of 
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interest (see below) the deflection angle is small, and propagating a distance I the proton is 
deflected by an angle (I / \) l / 2 \/ Rl. Thus, we may define an effective mean free path, the 
propagation distance over which large deflection occurs, by {IX) 1 / 2 / Rl = 1 and a diffusion 
coefficient for proton propagation, D = lc/3. For a propagation time t, protons are confined 
to a region of size d ~ (Dt) 1 / 2 = (ct/'5\) 1 l 2 R L . For t = t H « 10 10 yr, we have d ~ 
l(E/3 x 10 19 eVX^nGA^pJ -1 Gpc. The propagation distance is determined by the product 
.BA 1 / 2 . The upper limit on the intergalactic magnetic field implied by QSO Faraday rotation 
measurements, BX 1/2 < 1 nGMpc 1/2 ||13|,|1J], implies d > l(E/3 x 10 19 eV) Gpc. We conclude 
that the existence of a uniformly distributed inter- galactic magnetic field would have no effect 
on J max . 



C. Magnetic fields in large scale structures 

The discussion in the previous subsections shows that magnetic field in the source or 
a uniformly distributed inter-galactic field would not affect the upper bound / max - Could 
magnetic fields associated with large scale galaxy structures, i.e. clusters, filaments and 
sheets, affect the upper bound? 

Let us first consider galaxy clusters, where inter-galactic fields had been detected. The 
analysis of rotation measures of radio sources lying in the background of rich clusters implies 
the existence of strong fields between galaxies, B ~ lyuG with A ~ 10 kpc, in the central 



0.5 Mpc, cluster region [|Tj|. However, confinement of high energy protons in the cores of 
rich clusters would have little effect on our results, unless most of the high energy neutrino 
sources reside in the centers of rich clusters, which is not the case for either the hypothesized 
AGN or GRB neutrino sources. Moreover, high-energy protons can not be confined even 
in the central regions of a rich cluster, since for B ~ 1/xG and A ~ 10 kpc we have d ~ 
10(£/3 x 10 19 eV) Mpc over a Hubble time. 

We next turn to large-scale filaments/sheets. Kulsrud et al. has recently sug- 
gested that magnetic fields could be amplified by turbulence associated with the forma- 
tion of large scale filaments and sheets to near equipartition with turbulent kinetic energy. 
For characteristic turbulent velocities of Vt ~ 300 km/s on ~ 1 Mpc scale, and charac- 
teristic filament/sheet densities of n ~ 10~ 6 cm -3 , this scenario predicts magnetic fields 
B ~ 0.1(n/10 -6 cm _3 ) 1,/2 (v i /300 km/s) fiG in the high density large scale filaments/sheets, 
with coherence length L ~ 1 Mpc, comparable to the filament (sheet) diameter (thickness). 
It is not clear whether the suggested scenario for increase in magnetic field strength and 
coherence length to equipartition with the largest turbulent eddies can be realized. Further- 
more, even for a turbulent velocity of order several hundreds km per sec, the turn-around 
time of a ~ 1 Mpc eddy is longer than the Hubble time, and it is therefore not clear whether 
equipartition with the largest scale is achievable. Nevertheless, we consider this scenario 
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here since it is consistent with the upper limit, B < 1 /zG, implied for a field coherent over 
~ 1 Mpc inside high density large scale structures by QSO rotation measures [IT]. Although 
the Larmor radius of a ~ 10 19 eV proton is smaller than L ~ 1 Mpc, confinement of par- 
ticles would require a special field configuration. Even if such a configuration is produced 
by the random turbulent motions, which seems unlikely, variation of the field over scale L 
gives rise to particle drift velocity, Vd ~ Rlc/L, and therefore to the escape of particles on 
time scale t e ~ L 2 /cR L = 10 7 (L/1 Mpc) 2 (£/0.1//G)(.E/3 x 10 19 eV) yr. Since t e < t H , the 
hypothesized large-scale structure magnetic fields can not affect the bound on neutrino flux. 

Finally, we note that several authors have recently considered cosmic-ray proton propa- 
gation in an hypothesized large scale magnetic field, associated with our local super-cluster, 
of 0.1/zG strength and 10 Mpc coherence length corresponding to a hypothetical local tur- 
bulent eddy of comparable size |I8|JT9[| . For these parameters as well we have t e C and 
therefore even this field structure would not affect the neutrino bound. Nevertheless, two 
points should be made. First, if cosmic-rays are confined to our local super-cluster, one 
would expect the local cosmic-ray flux to be higher than average (since the production rate 
averaged over Hubble time should be higher than average in over-dense large scale regions), 
implying that the upper bound on the neutrino intensity is lower than J max derived here. 
Second, it is hard to understand how the hypothesized magnetic field structure could have 
been formed. The overdensity in the local super-cluster is not large and an equipartition 
magnetic field of strength 0.1/iG therefore corresponds to a turbulent velocity Vt ~ 10 3 km/s. 
A turbulent eddy of this velocity coherent over tens of Mpc is inconsistent with local pecu- 
liar velocity measurements (e.g. for a recent review). Moreover, the corresponding eddy 
turn around time is larger than the Hubble time. 



IV. AGN JET MODELS 

We consider in this section some popular models for neutrino production in which high- 
energy neutrinos are produced in the jets of active galactic nuclei Q. In these models, the 
flux of high-energy neutrinos received at Earth is produced by "blazars" , AGN jets nearly 
aligned with our line of sight. Since the predicted neutrino intensities for these models 
exceed by typically two orders of magnitude the upper bound, Eq. (fj), based on observed 
cosmic ray fluxes, it is important to verify that the models satisfy the assumption on which 
Eq. (g) is based, i.e. optical depth < 1 to p — 7 interaction. 

The neutrino spectrum and flux are derived in AGN jet models on the basis of the 
following key considerations. It is assumed that protons are Fermi accelerated in the jet to 
high energy, with energy spectrum dN p /dE p cx E~ 2 . For a photon spectrum dN^/dE^ oc 
E~ 2 , as typically observed, the number of photons with energy above the threshold for 
pion production is proportional to the proton energy E p (the threshold energy is inversely 
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proportional to E p ). This implies that the proton photo- meson optical depth is proportional 
to E p , and therefore, assuming that the optical depth is small, that the resulting neutrino 
spectrum is flatter than the proton spectrum, namely dN v /dE v oc E^ 1 , as shown in Fig. 
1. The spectrum extends to a neutrino energy which is ~ 5% of the maximum accelerated 
proton energy, which is typically 10 19 eV in the models discussed. 

The production of charged pions is accompanied by the production of neutral pions, 
whose decay leads to the emission of high-energy gamma- rays. It has been claimed pi 



that the observed blazar emission extending to ~ 10 TeV supports the hypothesis 
that the high-energy emission is due to neutral pion decay rather then to inverse Compton 
scattering by electrons. Thus, the normalization of the neutrino flux is determined by 
the assumption that neutral pion decay is the source of high-energy photon emission and 
that this emission from AGN jets produces the observed diffuse 7-ray background, $ 7 (> 
lOOMeV) = 10 _8 erg/cm 2 ssr [p3 |. Under these assumptions the total neutrino energy flux is 
similar to the 7-ray background flux (see Fig. 1). 

In the AGN jet models discussed above, the proton photo-meson optical depth r P7 at 
E p < 10 19 eV is smaller than unity. This is evident from the neutrino energy spectrum shown 
in Fig. 1, which is flatter than the assumed proton spectrum at E p < 10 19 , as explained 
above. In fact, it is easy to see that these models are constrained to have r p7 < 10~ 3 at 
E p ~ 10 19 eV. The threshold energy of photons for pair-production in interaction with a 
1 TeV photon is similar to the photon energy required for resonant meson production in 
interaction with a proton of energy E p = 0.2GeV 2 /(0.5MeV) 2 x 1 TeV = 10 18 eV. Emission 
of ~ 1 TeV photons from blazars is now well established P^fl , and there is evidence that 
the high-energy photon spectrum extends as a power-law at least to ~ 10 TeV |22J. This is 
the main argument used f2]| in support of the hypothesis that high-energy emission from 
blazars is due to pion decay rather than inverse Compton scattering. The observed high- 
energy emission implies that the pair-production optical depth for ~ 1 TeV photons is small, 
and that r P7 < 10 _4 (-E p /10 18 eV), since the cross section for pair production is ~ 10 4 times 
larger than the cross section for photo-meson production. This result guarantees that the 
upper bound J max on the neutrino intensity is valid for AGN jet models. 



V. GAMMA-RAY BURSTS 

In the GRB fireball model for high-energy neutrinos, the cosmic ray observations are 
naturally taken into account and the upper limit on high energy neutrino flux is automati- 
cally satisfied. In fact, it was the similarity between the energy density in cosmic ray sources 
implied by the cosmic ray flux observations and the GRB energy density in high energy pro- 
tons that led to the initial suggestion that GRBs are the the source of high energy protons. 
Just as for AGN jets, the GRB fireballs are optically think to 7- p interactions that produce 
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pions but-unlike the AGN jet models-the GRB model predicts a neutrino flux that satisfies 
the cosmic-ray upper bound discussed in Section II. 



A. Neutrinos at energies ~ 10 eV 



In the GRB fireball model Po| , which has recently gained support from GRB after- 
glow observations [2B|, the observed gamma rays are produced by synchrotron emission of 
high-energy electrons accelerated in internal shocks of an expanding relativistic wind, with 
characteristic Lorentz factor T ~ 300 In this scenario, observed gamma-ray flux vari- 



ability on time scale At is produced by internal collisions at radius m V cAt that arise 



from variability of the underlying source on the same time scale p8| . In the region where 
electrons are accelerated, protons are also expected to be shock accelerated, and their photo- 
meson interaction with observed burst photons will produce a burst of high-energy neutrinos 
accompanying the GRB ||. If GRBs are the sources of ultra-high-energy cosmic-rays [29,30 
then the expected GRB neutrino intensity is || 



1 

2' 

1.5 x 10~ 9 (J^J mmil^E^/E^GeVcm-h-hi- 1 , E h v « 10 14 eV. (6) 



Here, f n is the fraction of energy lost to pion production by high-energy protons. The 
derivation of f n and E h v and their dependence on GRB model parameters is given in the 
appendix [Eqs. ( |A2|) , (|A1|) 1. The intensity given by Eq. @ is ~ 5 times smaller than that 



given in Eq. (8) of ref. H, due to the fact that in Eq. (8) of ref. || we neglected the 
logarithmic correction In (100) = 5 of Eq. (0). 

The GRB neutrino intensity can be estimated directly from the observed gamma-ray 
fluence. The Burst and Transient Source Experiment (BATSE) measures the GRB fluence 
F 7 over a decade of photon energy, ~ O.lMeV to ~ lMeV, corresponding to half a decade of 
radiating electron energy (the electron synchrotron frequency is proportional to the square of 
the electron Lorentz factor). If electrons carry a fraction f e of the energy carried by protons, 
then the muon neutrino fluence of a single burst is E^dN u /dE u rs 0.25(/ 7r // e )F 7 /ln(3). 
The average neutrino flux per unit time and solid angle is obtained by multiplying the 
single burst fluence with the GRB rate per solid angle, ~ 10 3 bursts per year over Arc sr. 
Using the average burst fluence _F 7 m 6 x 10~ 6 erg/cm 2 , we obtain a muon neutrino intensity 
El$> v ~ 3 x 10 _9 (/ 7r // e )GeV/cm 2 ssr. Recent GRB afterglow observations typically imply 



f e ~ 0.1 p6 |, and therefore f w /f e ~ 1. Thus, the neutrino intensity estimated directly from 
the gamma-ray fluence agrees with the estimate @ based on the cosmic-ray production 
rate. 
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B. Neutrinos at high energy > 10 eV 



The neutrino spectrum (|6]) is modified at high energy, where neutrinos are produced by 
the decay of muons and pions whose life time r Mj7r exceeds the characteristic time for energy 
loss due to adiabatic expansion and synchrotron emission The synchrotron loss time 

is determined by the energy density of the magnetic field in the wind rest frame. For the 
characteristic parameters of a GRB wind, the muon energy for which the adiabatic energy 
loss time equals the muon life time, i?" is comparable to the energy E s at which the life 
time equals the synchrotron loss time, r*. For pions, E% > E%. This, and the fact that 
the adiabatic loss time is independent of energy and the synchrotron loss time is inversely 
proportional to energy, imply that synchrotron losses are the dominant effect suppressing 
the flux at high energy. The energy above which synchrotron losses suppress the neutrino 
flux is 



Here, L 7 = 10 51 L 7j 5ierg/s is the observed gamma-ray luminosity, At = lAt ms ms is the 
observed GRB variability time scale, E h ~ lMeV is the observed GRB photon break energy, 
T = 300T 300 , and £ e and £b are the fractions of GRB wind luminosity carried by electrons 
and magnetic fields. The observational constraints on these parameters are discussed in the 
appendix. At neutrino energy E v ^> E s v , the probability that a pion (muon) would decay 
before losing its energy is approximately given by the ratio of synchrotron cooling time to 
decay time t^/tv^) = [E„/E* , B v \]~ 2 , and the intensity of Eq. is suppressed by a 
similar factor, 



Since the wind duration, i.e. the time over which energy is released from the source, 
is T ~ 1 s, internal shocks may occur due to variability on time scale St larger than the 
source dynamical time, At ~ 1ms < St < T ~ 1 s. Collisions due to variability St > At 
are less efficient in producing neutrinos, f m oc St^ 1 , since the radiation energy density is 
lower at larger collision radii, leading to a smaller probability for photo-meson interaction. 
However, at larger radii synchrotron losses cut off the spectrum at higher energy, E s (St) oc St. 
Collisions at large radii therefore result in extension of the neutrino spectrum of Eq. (||) 
to higher energy, beyond the cutoff energy Eq. (|7|), and therefore yield E%$ u oc E^ 1 for 
E v > E^(At), since f n oc St^ 1 oc This extension is shown in Fig. 1. We note, 

that on time scale At ~ Is the expanding wind is expected to interact with the surrounding 
medium, driving a relativistic shock into the ambient gas. Protons are expected to be 
accelerated to high energy in this shock, the "external" shock, as well. The neutrino intensity 






1 -2 



GeVcm 2 s 1 sr 



E v ^>E t 
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and spectrum produced in the external shock are given by equations (|(||X2|J7]) with At ~ Is. 
Due to the low efficiency of the external shock, /„. ~ 10 -4 , its contribution to the neutrino 
flux is small. Note that as the external shock expands through a larger mass of the ambient 
gas it decelerates, and therefore on time scale At 3> Is the shock Lorentz factor is not large 
enough to allow acceleration of protons to high energy. 



C. Comparison with other authors 

In agreement with Rachen & Meszaros ||, we find that the neutrino flux from GRBs 
is small above 10 19 eV, and that a neutrino flux comparable to the 7-ray flux is expected 
only below ~ 10 17 eV. Our result is not in agreement, however, with that of ref. ||31|| , where 
a much higher flux at ~ 10 19 eV is obtained based on the equations of ref. 0, which are 
the same equations as used here. There is a numerical error in the calculations of ref. \\31\f \. 
Finally, we note that the highest energy to which protons can be accelerated increases with 
the collision radius £;™ ax oc dt 1 / 3 pOf| , and while _E™ ax ~ 10 20 eV for 5t ~ At ~ 1ms, collisions 



at larger radii, 6t ~ 0.03 s, are required to allow acceleration to the highest observed energy, 
~ 3 x 10 20 eV. In agreement with Rachen & Meszaros ||, we find that at this radius the 
neutrino spectrum produced through photo-meson interactions extends to ~ 10 18 eV [see 
Eqs. (0), ( Al|| . There is no contradiction, however, between production of high-energy 



protons above ~ 3 x 10 20 eV and a break in the neutrino spectrum at ~ 10 16 eV, since the 
efficiency of neutrino production at collision radius corresponding to 5t ~ 0.03 s is small and 
most of the flux is produced by collisions at smaller radii. 



VI. DISCUSSION 

We have shown that cosmic-ray observations set a model-independent upper bound of 
E 2 Q U < 2 x 10~ 8 GeV/cm 2 ssr to the intensity of high-energy neutrinos produced by photo- 
meson interaction in sources of size not much larger than the proton photo-meson mean-free- 
path, e.g. AGN jets and GRBs (see Fig. 1). This limit cannot be avoided by hypothesizing 
evolutionary effects of the sources (see discussion in §IIc) or by invoking magnetic field 
scenarios (see discussion in §111). 



1 The parameters chosen in |3T 



are L 7 = 10 50 erg/s, At = 10s, and T = 100. Using our 



equation (4) of ref. ||, which is the same as Eq. (|A2|) of the present paper, we obtain for 
these parameters f n = 1.6 x 10~ 4 , while the author of |TlJ obtains, using the same equation, 
U = 0-03. 
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Of possibly even greater interest to photon astronomers, we have shown that the cosmic 
ray measurements rule out the current version of theories in which the gamma-ray back- 
ground is due to photo-meson interactions in AGN jets. 

The neutrino flux predictions of AGN jet models are based on two key assumptions, 
namely that AGN jets produce the observed gamma-ray background and that high-energy 
photon emission from AGN jets is due to decay of neutral pions produced in photo-meson 
interactions of protons accelerated in the jet to high energy. Since the neutrino flux predicted 
by these assumptions is two orders of magnitude higher than the upper bound allowed by 
cosmic ray observations (see Fig. 1), at least one of the key assumptions is not valid, 
presumably the assumption that high-energy photon emission from AGN jets is due to 
photo-meson interactions. This conclusion is supported by multi-wavelength observations 
of the blazar Mkn 421, which show contemporaneous strong variability at TeV and X-ray 
energies with little evidence for GeV and optical variability p2|j33|] . This behavior suggests 
that the high-energy photon emission is due to inverse-Compton scattering by relativistic 



electrons 32 



Even if AGNs produce part of their high energy emission through 7 — p interactions, the 
upper bound derived here implies an upper bound of ~ 1 detected AGN neutrino per year 
in a high-energy neutrino telescope with an effective area of 1km 2 [|T[. The discussion in this 
paper shows that there is no observational motivation that would lead one to expect that 
the flux of high-energy neutrinos from AGNs will be measurable in a km 2 telescope. 

The upper bound Eq. (|2|) also applies to the intensity of high-energy neutrinos that 
may be produced through the decay of charged pions created by p — p — > n ± X (rather 
than p — 7) interactions, as long as the p — p optical depth in the source is not large. At 
present, predictions of high-energy neutrino flux based on such models are not available in 
the literature. 

The cosmic-ray flux below 3 x 10 18 eV is steeper than at higher energy. This is most 
likely due to a contribution to the cosmic-ray flux at lower energy from Galactic sources 
of heavy ions. This view has recently gained support from the detection in the Fly's Eye 
data of a small but statistically significant enhancement of the flux of cosmic-rays in the 
energy range of 2 x 10 17 eV to 3 x 10 18 eV along the Galactic plane ||. Extra-galactic sources 
of cosmic-rays may therefore exist that produce cosmic-rays with energy < 3 x 10 18 eV at a 
rate higher than given by Eq. ([[]). While we have no observational evidence for the existence 
of such sources, we can not rule them out based on cosmic-ray observations, and they may 
produce a flux of < 10 17 eV neutrinos which is higher than the upper limit implied by Eq. 
(0). Note that this argument does not affect the validity of the upper bound (0) for AGN 
models, since the neutrino emission from these sources peaks at ~ 10 18 eV. 

The neutrino flux predicted by the GRB model is consistent with the upper bound 
derived here. The intensity estimate we give here, Eq. (^), is ~ 5 times smaller than that 
we gave in ref . |3j , where the logarithmic correction of Eq. ([!]) was neglected. The intensity 



13 



calculated here implies a detection rate of ~ 20 neutrino induced muons per year for a 1 km 2 
detector (over 4tt sr). As discussed in ||, one may look for neutrino events in spatial and 
temporal coincidence (on a time scale of seconds) with GRBs. 

The GRB neutrino spectrum is consistent with the secondary-particle cooling constraints 
derived by Rachen & Meszaros 0. The neutrino flux above ~ 10 16 eV is suppressed, but 
this is also consistent with the acceleration of protons to > 3 x 10 20 eV (see Section Vc). 

Finally, we note that the GRB neutrino flux discussed here is the flux produced in situ, 
i.e. within the source. The energy loss of high-energy protons, > 5 x 10 19 eV, through 
photo-meson production in interaction with microwave background photons would lead to 
a background neutrino intensity (which will not be temporally associated with GRBs) com- 
parable to the upper bound shown in Fig. 1 at E v ~ 10 18 eV. This flux of high-energy 
neutrinos should exist regardless of the nature of the high-energy proton sources (assuming 
that these sources are indeed extragalactic, see f34|). 
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APPENDIX A: NEUTRINO PRODUCTION IN GRBS 



GRBs are possible sources of high-energy cosmic-rays p9| , |30| , which may account for the 
observed extra-Galactic high-energy proton flux 0,0. In the GRB fireball model ||25|| , which 



has recently gained support from GRB afterglow observations []26| , the observed gamma rays 
are produced by synchrotron emission of high-energy electrons accelerated in internal shocks 
of an expanding relativistic wind. The hardness of the observed spectrum, which extends to 
~ 100 MeV, requires wind Lorentz factors V ~ 300 [T27[]. In this scenario, observed gamma- 
ray flux variability on time scale At corresponds to internal collisions at radius »s T 2 cAt, 
which arise from variability of the underlying source on the same time scale [28|. Rapid 
variability time, 



1 ms, observed in some GRBs ||35|| , and the fact that a significant fraction 
of bursts detected by the Burst and Transient Source Experiment (BATSE) show variability 
on the smallest resolved time scale, ~ 10 ms |36|], imply that the sources are compact, with 
linear scale ro ~ 10 7 cm and characteristic dynamical time ~ 1 ms. 

In the region where electrons are accelerated, protons are also expected to be shock 
accelerated, and their photo-meson interaction with observed burst photons will produce 
a burst of high-energy neutrinos accompanying the GRB ||. The neutrino spectrum is 
determined in this model by the observed gamma-ray spectrum, which is well described by 
a broken power-law, dN^/dE^ oc Er: 13 with different values of (3 at low and high energy 
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37 1 . The observed break energy (where (3 changes) is typically E ~ lMeV, with /3 ~ 1 at 



energies below the break and (3 ~ 2 above the break. The interaction of protons accelerated 
to a power-law distribution, dN p /dE p oc E~ 2 , with GRB photons results in a broken power 
law neutrino spectrum §, dN v /dE u oc E~ p with /3 = 1 for E v < E%, and (3 = 2 for > ^ 
(see Fig. 1). The neutrino break energy E h v is fixed by the threshold energy of protons for 
photo-production in interaction with the dominant ~ 1 MeV photons in the GRB, 

El « 5 x 10 14 r2 00 (^/lMeV)- 1 eV, (Al) 

where Y = 300r 30 o- 

The normalization of the flux is determined by the efficiency of pion production. As 
shown in |J , the fraction of energy lost to pion production by protons producing the neutrino 
flux above the break, is essentially independent of energy and is given by 

U = 0.20-^- . (A2) 

j7T (E*/lMeV)Tj 00 At ms 1 J 

Here At = lAt ms ms and L 7 = 10 51 L 7i 5ierg/s is the observed gamma-ray luminosity. The 
values of T and At in Eq. (|A2]) are determined by the hardness of the 7-ray spectrum 



and by the flux variability. These parameters are also constrained by the fact that the 
characteristic observed photon energy is ~ 1 MeV. Internal collisions are expected to be 



"mildly" relativistic in the fireball rest frame pq| , i.e. characterized by Lorentz factor 
7i — 1 ~ 1, since adjacent shells within the wind are expected to expand with similar 
Lorentz factors. The internal shocks would therefore heat the protons to random velocities 
(in the wind frame) 7 P — 1 ~ 1. The characteristic frequency of synchrotron emission is 
determined by the characteristic energy of the electrons and by the strength of the magnetic 
field. These are determined by assuming that the fraction of energy carried by electrons is 
£ e , implying a characteristic rest frame electron Lorentz factor 7 e = £ e (m p /m e ), and that a 
fraction of the energy is carried by the magnetic field, implying Aur^cT 2 B 2 /8n = ^bL 
where L is the total wind luminosity. Since the electron synchrotron cooling time is short 
compared to the wind expansion time, electrons lose their energy radiatively and L w £ 7 /£ e - 
The characteristic observed energy of synchrotron photons, E h = Th^eB /m e c, is therefore 



L 1 ' 2 



E b , - ^^^Ag-MeV. (A3) 

At present, there is no theory that allows the determination of the values of the equipartition 
fractions £ e and However, for values close to equipartition, the model photon break 
energy E^ is consistent with the observed E^ for T = 300 and At = 1 ms. 
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